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TECHNICAL NOTE NO. 1014 

STRESSES IN AND GENERAL INSTABILITY OE MONOCOQUE CYLINDERS 

WITH CUTOUTS 

II - CALCULATION OP THE STRESSES IN A CYLINDER 
WITH A SYMMETRIC CUTOUT 

By N. J. Hoff, Bruafi A. Boley, and Bertram Klein 

SUMMARY 


A numerical procedure is presented for the calculation 
of the stresses in a monocoque cylinder with a outout. In 
the procedure the structure is broken up into a great many 
units; the forces in these units corresponding to specified 
distortions of the units are calculated; a set of linear 
equations is established expressing the equilibrium con- 
ditions of the units In the distorted state; and the 
simultaneous linear equations are solved. A fully worked 
out numerical example, corresponding to the application of 
a pure bending moment, gave results in good agreement with 
experiments carried out earlier at the Polytechnic Institut 
of Brooklyn. 


INTRODUCTION 

Actual airplanes differ greatly from the idealized 
structures that underlie most theoretical analyses. The 
reason for these deviations can be found in the great 
difficulties involved in applying the theory of elasticity 
to the Irregular and complex structural parts of airplanes. 
It is believed that the most promising approach to these 
complex problems is the one in which the structure is 
imagined to be broken up into a great number of "unita," 
the forces in these units corresponding to specified dis- 
tortions of the units are calculated, a set of linear 
equations is established expressing the equilibrium con- 
ditions of the units in the distorted state, and the 
simultaneous linear equations are solved. 
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The set of linear equations, excluding the load terms, 
forms what Southwell (reference l) called the "operations 
table." In Southwell's relaxation procedure the equations 
are solved by a method of step— by— step appr oxima.t ions . In 
the past 2 years a considerable amount of work has been done 
at the Polytechnic Institute of Brooklyn in applying 
Southwell's method to the stress analysis of reinforced thin- 
walled structures. It was easily possible to establish a 
rapidly converging procedure in the case of stiffened panels 
(references 2 and 3;. In the case of ring problems (refer- 
ences 4 and 5) the convergence was found to be poor as a 
rule, and suggestions were made for solving the equations 
either directly by matrix methods, or by a procedure denoted 
as the "growing unit" method. 

The problem of the calculation of the stresses in a 
reinforced monocoque cylinder combines the two elements 
discussed in the earlier reports, namely, reinforced panels 
and rings. The authors were unable to devise a pure step- 
by-step procedure that would lead to a solution of the 
equations represented by the operations table with a reason- 
able expenditure of work and time. On the other hand, the 
solution can be found comparatively easily if the operations 
table is set up with the aid of the expressions developed 
in this report and the equations which the operations table 
represents are solved by matrix methods. A fully worked 
out numerical example gave results in good agreement with 
the tests described in reference 6. 

This investigation, conducted at the Polytechnic 
Institute of Brooklyn, was sponsored by and conducted with 
the financial assistance of the National Advisory Committee 
for Aeronautics. 


SYMBOLS 


a distance between rings measured in the x— 

direct ion 

& * * c , d , e , f ,g, 
h,j ,k , 1 ,m , n ,p , 
q,r ,s ,t ,u,r' , 

s' , t' portions of operations table 

a n«^n» c n Fourier coefficients 

A cross-sectional area of a stringer plus its 

effective width of sheet 
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4 

4 




H,J,E,0,P point* at intersection of rings and stringers 

A,B,C,D,E,3r.G, 

H.J.K.l.h.N.P, 

Q.R.S.T.U.Y.V,,, 
rf,X,Y,Z,U« ,Y« , 

Vi.rf' ,Y‘ ,Z‘ portions of operations table 

A-A, B-B , 0-0, 

D-D, A’— A’ , 

B * — B 1 , G 1 — C * , 

D D ' rings 


AB, BO, OD, 

ED 1 field* between rings 

E Young* s modulus 

£ g shear stress 

0 shear modulus 


1r moment of inertia of ring cross section plus 

its effective width for bending in its 
plane 

L distance between stringers measured along 

the circumference, that is, developed 
length of ring segment between adjacent 
stringers 

M bending moment acting in a transverse section 

of the cylinder 

n index 


bending moment acting in the plane of a ring 


PiiPa 


unknown coefficients 


P 

<1 


resultant force acting in a transverse section 
of the cylinder 

shear flow 


r 


radius of cylinder 



NA.CA. TN Ho. 1014 


4 


B. radial fora© 

SS line of syaoetry 

t sheet thickness 

T tangential force 

u tangential displacement of an intersection point 

of a ring and stringer 

v radial displacement of an intersection point of 

a ring and stringer 

w rotation of a section of ' a ring in its own plane 

X axial force 

a n ,a r ,a t coefficients used in the calculation of forces 

and moments caused by the shear flow existing 
in a panel 

0 angle subtended by ring segment 

V s ect ion— length parameter 

T = G-t/2 

T| vertical downward translation of a ring as a rigid 

body 

A = GtL/a 

1 axial displacement; also the ratio of effective 

shear area of a section to the actual area 
of the section 

cp angular coordinate 

rigid body rotation of a ring 

Q - G-t a/4L 

The symbols used to denote influence coefficients are 
defined in the following manner: 

( ab ) stands for the force or moment a caused by a unit 
movement in the direction of b (which direction is that of 
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the force B or T, or of the moment N). Thus (nil) is the 
moment due to a unit rotation; while (^r) is the tangential 
force arising from a unit radial displacement. Further, to 
distinguish the reactions at the fixed end from those at the 
movable end. the subscripts F and M are employed. Conse- 
quently, (nS)ji is the moment arising at the fixed end of 

the curved bar as a result of a unit tangential displace- 
ment of the movable end; while stands for the tan- 

gential force at the movable end due to a unit tangential 
displacement of that end. 

It should be noted that all the reactions considered 
here are acting from the support upon the curved bar. 


CALCULATION OF. THE DEFORMATIONS AND STRESSES IN THE 

SIMPLIFIED CYLINDER 
Assumptions Regarding the Structure 


The actual monocoque cylinder (fig. la) contains 1 6 
stringers and 8 rings, including the 2 end rings. The cut- 
out In the cylinder extends over three ring fields and two 
stringer fields. There are, therefore, 106 panels in the 
structure not counting the 6 panels cut out. Since the 
"unit" of the structure is the panel and its bordering 
stringer and ring segments. It appears that there are too 
many units to permit a calculation of the displacement and 
stresses with a reasonable amount of work. Hence it was 
decided that the actual structure should be replaced by the 
simplified one shown in figure lb. 

It was anticipated that the effect of the cutout upon 
the stress distribution would be important only in the 
neighborhood of the cutout. This is the justification for 
choosing small size units close to and large size units 
farther away from the cutout. In order to avoid a change 
in the total amount of stringer cross section in the 
cylinder, the sizes of the stringers bordering large panels 
were increased, as may be Been from the data contained in 
figure lb. The total width of sheet between adjacent 
stringers was considered effective in carrying normal stress 
in the axial direction and was distributed evenly between 
the adjacent stringers. Because of this assumption the 
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results of the calculations can "be anticipated to agree 
with experiment only when the loads applied are so small 
that the panels of sheet are in an unbuckled state. When 
calculations are to he carried out for higher loads, a 
smaller value should he assumed for the effective width. 

The effective width of sheet acting with the ring was 
taken equal to the width of the ring. 

Because of the symmetry of both structure and loading 
only one— quarter of the monocoque cylinder was considered 
in the calculations. 

In regard to the mechanical properties of the elements 
of the structure the following assumptions were made: 

1. The stringers have only extensional rigidity; they 
are very weak in bending and torsion. 

3. The rings are resistant to shearing and bending in 
their plane of curvature and to extension, but are very 
weak in bending out of their plane and in torsion. 

3. The sheet is resists.nt only to shearing deformations. 
The extensional rigidity of the sheet is taken into account 
through the assumption of an effective width the area of 
which is added to the cross-sectional area of stringers 
and rings. A conseauence of this assumption is that the 
shear stress must be constant over any single panel. 

These assumptions are believed to represent the essential 
features of the elements of the structure. The agreement 
between the results of the calculations and the experiments 
substantiates this belief. 


The Unit Problem 

1. Axial displacement t at point A .— When point A 

in figure 2 is displaced a distance | axially (in the x— 
direction), the force g ^ r exerted by the stringer segment 

AD upon the constraint at D is 




(la) 


Xp str 


= ( s ^str/ a H 
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where A g ^. r is the cross-sectional area of the 
plus its effective width of sheet. At the same 
constant shear flow q^| will he caused in the 


8 tr inger 
t ime a 
panel : 


q A| = ( Gt / 2L ^ (lb) 

In figure 2 the shear flow is indicated as acting from the 
sheet upon the edge reinforcements. As in reference 2 one- 
half the shear force transmitted from the sheet to the edge 
stringer is considered to he acting upon each of the con- 
straints at the end points of the stringer. Consequently, 
the axial force Xjj transmitted hy shear to the con- 

straint at D is 


X D sh= -(<Ha/4L)£ 


( lc) 


The total axial force X-n acting upon the constraint at 
D is 


X D = C(BA str /a) - (Gta/4L)]| (l) 


The shear flow q^_ acting upon the ring segment AB 
causes forces and moments to act upon the constraints at 
A and B which can he computed from the data given in 
reference 5. It should he remembered, however, that in 
this reference ths quantities listed are forces and moments 
acting from the constraints upon the ring segment at the 
end point toward which the shear flow is directed, and the 
sighs correspond to the beam convention (fig. ?), At the 
other end the same numerical values apply, hut the signs are 
inverted. Since in figure 2 the shear flow is directed to- 
ward point A, and since the forces and moments acting upon 
the constraint are sought, a multiplication hy —1 must he 
carried out. The forces and moments acting upon the con— 
straint at A are, therefore, according to the beam con- 
vent ion 


T A * - a t Lt U “ ~ 

R A = ” a r Lt lA = ~ a r (Gt/2)S 


> ( 3a) 


N 


A 


“ a n L °-A 


_L(Gt/3)| 


= — a. 
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Similarly at point B 


H]3 = cc, r Lq A = a r (Gt/2)| 

N b = a n L*q A = a a L(&t/2)! 


"A 


> 


J 


( 2b) 


If these equations are rewritten to agree with the 
frame convention (fig. 3), the following is obtained: 


T k = cx t (&t/2)! 
E a = -a p (&t/2)| 
N a - a£(St/2)| 

T 3 * « t (Gt/2)| 
B b = « r (Gt/2)e 

K B = a n LC&t/2)| 


■'l 



is 


Since the shear 
equal and opposed 


flow transmitted to ring segment CD 
to that transmitted to ring segment AB, 



It is easily verified' on the basis of the principles 
developed in reference 2 that the axial forces acting upon 
the constraints at points A, B, and C are: 


X A = -[(EA gtr /a) + (Gta/4L)3| 
X B = (Gta/4 1)1 
X c = (Gt»/4L)t 


(5) 

( 6 ) 
(7) 
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2. Tangential Si solacement at Point A .- When point A 
is displaced tangentially through a distance u (in the pos- 
itive direction of the "beam convention) the forces and moments 
acting upon the ring segment at A and 3 can he computed 
from the data presented in reference 5. 


27 A = tt U yi 

R A “ 

= £t K u 

* ^ y u 

= nt-pU 







( Sa) 


(8h) 


At point A a positive u displacement according to 
the frame convention is equal to a negative u displacement 
according to the "beam convention. If a positive u dis- 
placement is now assumed according to the frame convention, 
the forces and moments are expressed as they are acting 
upon the constraint, and the signs are written in agree- 
ment with the frame convention, equations (8a) and (8b) 
become 


*A 


-t t K tt 




R k ~ rt M u 
N k = ~^M U 


> 


( 8 ) 


T 3 = ttj,u 


E B = rt ]T u 




(9 ) 


n b = at r u 
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At the same time the displacement u gives rise to a 
uniform shear flow in the panel. The direction of the shear 
flow acting upon the edge reinforcements of the panel is 
indicated in figure 4. Its magnitude can "be calculated with 
the aid of Maxwell’s reciprocal theorem. When A was dis- 
placed in the x— direction, the induced shear flow caused a 
tangential force at A of a magnitude 


®A * 

according to the first of equations (2). Hence when A is 
displaced tangentially a distance u, according to the 
reciprocal theorem the axial force caused at A must he 
equal to 


X x = a^( G-t /2 )u 


( 10 ) 


This axial force, however, can occur only if the 
magnitude of the shear flow <1 a.u ^7 


^Au 


a t ( Gt /a)u 


( 11a) 


The direction of the shear flow must be determined from 
equation (10). Since a^. is negative within the range of 

angles represented in the graphs of reference 5, the force 
exerted upon the constraint at A is directed toward the left 
in figure 4, Cons equently, the shear flow acting upon the 
edge reinforcements of the panel must be as shown in figure 
4. 


It may be seen that the direction of the shear flow is 
the same in figures 2 and 4. Equations (2) and (3) can be 
easily transformed to correspond to the values caused by 
the tangential displacement. 


T^ =s —a*( GtL/a)u 

R a « a r a t ( GtL/a)u 
s 

N a = — c^a^C&tl /a)u 


( 11 ) 
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= — a.®( GtL/a)u 

Eg = — <x r a^.( GtL /a )u / 

N b = -a^.o^ GtL S /a)u 

Equations (4) are valid here also 

The total forces and moments at A are obtained b 
ing up corresponding values in equations (s) and (ll) 

2 — \ 

T^ = — + a^.( GtL/a ) ]u 

a A = + cXpO^C GtL/a) 3u > 

*b — [ntg + ct, a a^.(GtL /a)3u 

Similarly the total forces and moments at B can 
calculated from equations (9) and (12): 

T b = [ttj, — a^( GtL/a) ]u 

R b = Crtj, - GtL/a) ]u 

R B = [nt y — a^a^C GtL S /a) 3 

At G the forces and moments are: 

jg 

T q = a^.( GtL/a)u 

R_ = a <x, ( G tL/a)u 
0 r t 

H- = a a. (GtL s /a)u 
u n t 

At D the forces and moments are: 

T d = a*(GtL/a)u 

= — a r a^. ( GtL/a)u 
= a n a t (GtX< /a)u 


> 

J 


> 


> 

- 


( 12 ) 


add- 


(13) 


be 


( 14) 


( 15) 


( 16 ) 
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The axial forces can "be calculated from the following equa- 
tions : 

X B = X 0 = -X D = -X A = — a t (Gt/2)u (17) 


3. Radial Displacement v and Rotation w at Point A. 

The calculations needed for determining the forces and mo- 
ments corresponding to these distortions are quite similar 
to those given under (2). The results of the calculations 
are presented in the diagrams to be discussed under (4). 

4. The four— panel problem, — In the general case any 

point belongs simultaneously to four different panels. A 
displacement of the point, therefore, will cause forces and 
moments to appear in four panels. These forces and moments 
can be calculated without difficulty from the results of the 
single— panel unit problem. For the convenience of the stress 
analyst the four— panel problem has been worked out and the 
results of the calculations are presented in figures 7 to 10. 
In figure 10a the sign convention is shown. 


When the sheet is in a buckled state in any particular 
panel, a reduced value should be used for G. If one or two 
of the four panels adjacent to a point are cut out, G should 
be put equal to zero for those panels. 


The Operations Table 

The operations table contains the forces and moments 
acting upon the constraints caused by the individual unit 
displacements. The individual items in the table are calcu- 
lated according to the principles presented in the preceding 
section. Each number in the operations table represents the 
value of the force quantity indicated at the left end of the 
row in which the number is located, caused by the unit dis- 
placement indicated at the top of the column in which the 
number is located. This arrangement differs from the one 
used in references 2 to 4 insofar as the headings of the 
columns in this table are those of the rows in the references 
mentioned, and vice versa. The sets of the headings of the 
columns and rows, however, are interchangeable since the 
operations table is symmetric with respect to its principal 
diagonal. The number of individual operations — that is, 
the number of degrees of freedom of the structure — can "be 
determined with the aid of the following considerations. 
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A s was stated earlier, it suffices to consider only- 
one— quarter of the entire structure in the calculations be- 
cause of the symmetry of both structure and loading. 

Points E, 1, M, 3?, and 0 of the end ring are considered 
rigidly fixed. In the tests described in reference 6 they 
were attached to a rigid frame. Points B, 0, D, P, G, and 
H are free to more axially, tangentially, and radially, and 
the sections of the rings at these points are free to rotate 
in the plane of the ring. Points A, 35, and J are free to 
move axially and radially, their other two types of motion 
being excluded because they are antisymmetric with respect 
to the plane of symmetry of the cylinder passing through 
points K, A, 0, E, and JT. 

Altogether, the system has 30 degrees of freedom. The 
displacements corresponding to them are arranged according 
to the following scheme in the operations table (table l). 
Pirst, all the axial displacements are listed, nine in 
number. They are followed by all the other displacements 
of each of the points of ring ABODE, arranged in the order 
of tangential displacement, radial displacement, and rotation. 
Altogether, there are 11 such operations if the antisymmetric 
distortions at points A and E are excluded, Pinally, the 10 
individual operations in the plane of ring PGHJ are listed. 

In the appendix it is shown by means of typical examples 
how the entries in the operations table are calculated. 


Calculation of Displacements 

As far as the loading is concerned, the forces acting 
in the end sections upon the individual stringers are not 
stipulated, but it is required that they add up to a pure 
bending moment acting in the vertical axial plane of symmetry 
of the cylinder. On the other hand, it is known that ring 
PG-HJ must remain plane during the distortions. The calcu- 
lation of the distortions of the structure is, therefore, 
carried out in the following manner. 

Plane PGHJ is assumed to be rotated about the horizontal 
axis in its plane passing through point H. The angle of 
rotation is defined by the assumption that the axial dis- 
placement of point J is 0.001 inch. 
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It was shown In reference 4 that the operations table 
represents a set of simultaneous linear equations. Tor 
instance, the first row in the operations table (table l) 
may be written in the form 

-113.9! + 38.9! b + 0.934v a + 11.9u B - 0.934v B + 0.305w B = 0 

As a second example, row 5 reads: 

1?I D “ 7196 E + 7.3| h + 394. 76 j 

— 11.7UJJ — 3.78t h — 5.05w H + 3.78Vj « 0 

In this equation the values of ! H and are 0 and 1, 

respectively, according to the assumptions made regarding the 
rotation of ring FGHJ . (in the operations table the unit 
axial displacement is 0.001 in.} Conseauently , the entry 
in column 9 and row 5 is a known quantity having the value 
394.7 X 1 « 394.7 pounds. It can be taken over to the right- 
hand side of the equation. Since the entry in column 8 and 
row 5 is zero because * 0, the equation corresponding 

to row 5 may be written as 

17| D - 719| b - 11.7u h - - 3 -78v h - 5.05v H + 3.78VJ = -394.7 

It is easy to see that the complete set of equations 
corresponding to the assumed rotation of ring FGHJ can be 
obtained by multiplying the figures listed in each column 
corresponding to an assumed axial displacement by the assumed 
value of the displacement, transferring the numbers obtained 
to the right-hand side of the equations, and adding up 
algebraically the numbers on the right-hand side of each 
equat ion . 

The terms contained in rows 6 to 9 of the operations 
table add up to the axial forces acting upon points F, G, H, 
and J , respectively. They need not be equated to sero since, 
because of the symmetry, equal and opposite forces originating 
from the omitted other half of the cylinder automatically 
balance them. For this reason, rows 6 to 9 must be omitted 
from the set of equations to be solved. They will be used 
later for establishing the nature of the external loading 
of the cylinder. 
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The set of the remaining 26 equations was solved ‘by- 
Doolittle’s method (references 7 and 8). In other words the 
26 displacements and rotations were calculated that correspond 
to zero resultant force and moment at each point of the struc- 
ture, It should he noted, however, that equations of equi- 
librium at points of the fixed ring KLMNO were not taken in- 
to account since the rigid fixation is capable of providing 
any reactions that are needed for equilibrium. Similarly, 
the equations of equilibrium of the axial forces at the 
middle ring FG-HJ were not considered because they were 
balanced by equal and opposite forces arising from the sym- 
metric other half of the cylinder, as mentioned before. The 
forces and moments arising in the two symmetric halves of 
the cylinder in the T, E, and N directions do not balance one 
another at the middle ring i'&HJ but add up, since they are 
equal In both magnitude and Bense, The conditions of equi- 
librium of these forces and moments are consequently included 
in the set of equations. 


The displacement quantities obtained follow; 


&A 


-0.1055 

^B 

3 

-0.4676 


ae 

-0,4105 

^D 

= 

0,0001 


3 

0,5735 


ss 

-0.9239 

^G- 

sc 

-0.7071 


= 

0 


= 

1.0000 

V A 

3 

2.0291 







U B 

= 

0.3233 

V B 

3 

-0.9576 

W B 

= 

-0.7573 

u c 

= 

-0.2924 

v c 

3 

-1.2685 

w c 

a t 

0,4316 

U D 

3 

-0.1146 

V D 

3 

0.5806 

W D 

sx 

-0.2831 


3 

1.3471 








3 

-0.3165 


3 

0.1872 

W P 

s 

-0.1772 

U G 

3 

-0.4035 

V G 

3 

-0.8034 

W Cr 

* 

-0.4348 

U H 

3 

-1.7210 

V H 

3 

-0.7722 

W H 

ae 

0.5665 

V 

J 

3 

-0.6707 








(18) 
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It may be noted that the unit displacement is 0.001 inch, 
and the unit rotation 0.001 radian. 

Substitution of the values obtained into the equations 
corresponding to the sixth to ninth rows of the operations 
table gives the resultant forces acting upon the constraints 
at points J to J in the axial direction. The forces in the 
stringers at these points can be obtained by multiplying by 
— 1 the forces calculated: 


Xy = —79.9084 pounds 
Xjj as —0,106 4 pound 


Xq = —89.1049 pounds 

( 19 ) 

Xj = + 170.9411 pounds 


The forces correspond to a bending moment 

M = 3077.7 inch-pounds (SOa) 

and a tensile force 

P a 1.8214 pounds (20b) 

Hence the rotation of the plane of ring FG-HJ undertaken 
corresponds to the application of a considerable bending 
moment and a very small tensile force. The tensile force 
can be eliminated by a suitable axial translation of the 
plane of ring FGHJ , 

The second part of the calculations consisted, therefore, 
of the determination of the distortions of the cylinder corre- 
sponding to an axial displacement of ring FG-HJ amounting to 
—0,001 inch. The simultaneous linear equations were set up 
in the same manner as before, except that now the values 


= ^ Q. = ^ h “ = ~ 1 


were used. This new system of equations was solved by 
Doolittle’s method. Those acquainted with the method will 
realize that this second solution involves comparatively 
little work if use is made of the solution of the first set. 
The results are: 
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£ A 

= 

-0.0993 


— 

-0.5039 


8 

-0.5803 


= 

-0.5728 

^23 

= 

-0.5725 

Ip 

3 

-1.0000 

^ a 

s 

1.0000 

t E 

= 

-1.0000 


= 

-1.0000 

V A 

= 

2.4041 







U B 

= 

0.4899 

V B 

s 

-0.4650 

"B 

3 

-0.7013 

U C 

s 

0.0674 

T o 

s 

-0.8912 

w 0 

3 

0.3489 

U D 

s 

0.0611 

V D 

3 

0.1735 

W D 

8 

-0.09 00 


8 

0.2400 







u r 

ss 

0 . 4055 


= 

0.3407 

W S> 

= 

-0.2371 

u o 

ss 

0.3176 

V 

Gr 

8 

0.7220 

w o 

3 

- 0.1459 

U H 

= 

-0.2376 

V H 

8 

-0.0716 

W H 

3 

0.2123 


= 

-0.4981 








> ( 21 ) 


J 


The unit displacement is again 0.001 inch, and the unit 
rotation 0.001 radian. 

Substitution of the values obtained into the equations 
corresponding to the sixth to ninth rows of the operations 
table gives the resultant forces acting upon the constraints 
at points S’ to J in the axial direction. The forces in the 
stringers at these points can be obtained by multiplying by 
— 1 the forces calculated: 


Xp = —87.6409 pounds 
Xg = —257.4269 pounds 


X_ = —126.0882 pounds 


X = —171.9484 pounds 
0 


( 22 ) 
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These forces correspond' to a resultant force 

P a —643,1044 pounds (23a) 

and a "bending moment 

M ** —1.821 X 10 inch— pounds ( 23b) 


Obviously, the two solutions — namely, those corre- 
sponding to the purs rotation and the pure translation, 
respectively, of ring PGHJ — can he combined in such a 
manner as to represent the two loading cases of pure bend- 
ing and pure compression. The solution may be obtained by 
solving in each case two simultaneous equations. Because 
of the great difference in the numerical values, however, 
it is quicker and just as accurate for practical purposes 
to correct for the effect of tension in the following 
manner in the loading case corresponding to bending: 

The tensile force caused by the rotation is 1,8214 pounds. 

The compressive force due to the translation is 643.1044 
pounds , 

A combination of the rotation undertaken with a trans- 
lation of 1.8214/643.1044 = 0.00283 units eliminates the 
tensile force and introduces an additive bending moment of 
less than one— hundredth of a percent of the original moment. 

The final pattern of distortions is sh».wn in figures 
11 to 14. 

In a similar manner, the loading case eorr esponding to 
pure compression can be dealt with. The dieplacements calcu- 
lated for pure translation must be combined with those calcu- 
lated for the pure rotation multiplied by the factor 
18.21/3077.7 = 0.00593. 

The final pattern of distortions corresponding to pure 
compression is shown in figures 15 to 18. 
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Calculation of the Stresses 

The average normal stress In a segment of a stringer 
between two adjacent rings is equal to the difference 
between the axial displacements of the end points of the 
segment times Young's modulus of the material divided by 
the original length of the segment. The average stress 
was calculated from the displacement values obtained as 
shown in the preceding section. The values of the stress 
are plotted against the distance of the stringer from the 
horizontal axis of the cylinder in figures 19 to 22. The 
curves shown correspond to either a pure bending moment of 
35,000 inch-pounds or to a pure compression of 1J386 pounds. 
The former value was chosen in order to permit a comparison 
with experimental results. 

Theory and experiment agree in obtaining a practically 
linear stress distribution over the major portion of both 
the complete and the cut sections. The straight lines, 
however, do not coincide. The reason for the discrepancy 
is the difference in the location of the centroids of the 
actual and the simplified monocoque cylinders. Deviations 
from the straight line occur in the neighborhood of the 
cutout. The nature and the magnitude of these deviations 
are practically the same in experiment and calculation. 

The shear stress in a panel depends upon the displace- 
ments and rotations occuring at all the four corners of the 
panel. If in figure 23 point A is displaced axially through 
a distance | A , the shear stress induced in the panel is 

f g * ~( 1/2 )(G/L)! a ( 24a) 

provided the positive sense of the shearing stress acting 
in the sheet is as shown in figure 23, It follows from 
equation ( 11a) and the sign convention adopted that the 
shear stress caused by a tangential displacement of an 
amount u A of point A is 

f s ** (24b) 

Similarly the shear stress caused by a unit radial displace- 
ment v A of point A is 

f a = a v (S/a) TjL 


( 24c) 
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Finally a rotation w A of the ring section at A gives rise 
to a shearing stress 


f s = a^,(G/a)v A 


(24d) 


Displacements and rotations at the other comers of 
the panel contribute similar quantities to the shear stress, 
hut care must he taken to use the proper sign. The total 
shear stress is 


^s = (d/2) (G-/L) (- + Sb + §£ " Sp)* ct t( < j/ a ) ( u A + ' a B~ u C" u D^ 


+ 0 ^.( 0 ^) ( -v A +v B -v c +v D )+a n L(G/a) (w^Vg-WQ-v^) (24) 


Substitution of the values of the constants and the 
displacement quantities yields the shear stress in any panel. 
Figures 24 and 25 contain shear stress distribution curves 
for pure bending and pure compression, respectively. A 
comparison of the calculated values with experimental ones 
is not well possible because of the limited number of 
measurements and because of the simplifying assumptions of 
the calculations. Experimental strain data are not avail- 
able for the full section. In the cut section the average 
measured value in the panels adjacent to the cutout was 
454 psi (fig. 36 of reference 6), while the calculated 
value was 4o.psl. This latter, however, was obtained as 
a small difference of large quantities and is not reliable 
for this reason. Measurements in the full section were 
made in a cylinder having the large cutout. The results 
(see fig. 40 of reference 6) indicate considerably higher 
stresses in the panels of the full section than in those 
of the cut section near the cutout. Moreover, a change 
in the sign of the shear also was observed. 

The calculation of the bending moments, shear forces, 
and tensile forces In the rings can be carried out accord- 
ing to the principles stated in reference 4, In the present 
case the maximum bending moment in the full ring at the edge 
of the cutout was found to be 2.64 inch -pounds when the load- 
ing of the cylinder was a bending moment of 35 ,,000 inch-pounds. 
The moment diagram is shown in figure 26. 
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The 2 . 64— inch— pound 'bending moment ia very insig- 
nificant as compared to the applied bending moment of 35,000 
inch— pounds. Nevertheless, it causes high stresses beeause 
of the small moment of inertia of the ring section. The 
maximum stress is 2,64 X 850 = 2240 psi according to the 
Mc/I formula, 

When the applied loading is a compressive force rather 
than a bending moment, the moments in the ring are found 
to be similar to those Just discussed. For this reason the 
moment diagram is not shown. 


DEVELOPMENT OP A STEP-BY-STEP APPROXIMATION PROCEDURE 

Basic Considerations 


The method of calculation of the stress distribution 
in a monocoque with a cutout presented in the first part of 
this report gave satisfactory results with a reasonable 
expenditure of work. It is felt, however, that with more 
complex structures — for instanoe, monocoque* with non— 
symmetric cutouts or monocoques with several cutouts — the 
operations table would become so large that its solution by 
the methods of matrix calculus might entail too much nu- 
merical work to be practicable. An effort is made in this 
part, therefore, to develop a procedure of step— by— step 
approximations suitable to cope with these complex problems. 

In order to simplify the presentation, the structure 
shown in figure 27 Is used in place of the actual monocoque 
cylinder. This structure will be referred to as the small 
cylinder. The vertical transverse plane of symmetry S— S is 
considered here as fixed in space and the points of the 
cylinder on the four rings are mqved relative to this fixed 
section. However, end rings A— A and A 1 — A* are rigid and can 
only undergo rigid body displacements. Because of the 
double symmetry of structure and loading it suffices to 
list only the forces and moments caused at points contained 
in one— quarter of the cylinder. 

In the actual calculations all four rings were first 
rotated as rigid bodies about their horizontal diameters. 

The amount and sign of rotation were defined by the stipu- 
lation that the intersection point of ring A— A with stringer 
1 be displaced a distance of 0.003 inch, that of ring B— B 
with stringer 1 a distance of 0,001 inch, both in the 
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negative x— direction. The axis of rotation was the horizontal 
diameter of the ring. Rings A 1 — A. 1 and B*— B* were rotated 
symmetrically. This distortion pattern corresponds to that 
prevailing in the complete cylinder (without the cutout) 
under the action of a pure "bending moment, provided rings 
A— A and A'— A' also undergo a rigid body translation verti- 
cally downward. The necessary amount of translation was 
determined from the requirement that at each point along 
ring B— B the forces and moments caused by the displacements 
had to add up to zero force and moment resultants. In 
practice, only one component force or moment had to be 
balanced out at any single point, after which all the other 
points were found to be automatically in equilibrium. 


As the next step the unbalanced forces and moments 
were calculated that arose in the structure of figure 27 
when the displacements determined in the preceding para- 
graph for the complete cylinder were applied to the cut cyl- 
inder of figure 27. Because the loading did not involve 
shear forces, unbalances were found only in the x— direction 
and at points along the edges of the cutout. 

It was anticipated that the distortion pattern would be 
influenced materially by the cutout only in the neighborhood 
of the cutout. Consequently, additional displacements would 
have to be undertaken only at the intersection points of 
stringers 7, 8, and 9 with ring B— B , This restriction 
materially decreased the amount of work involved in the 
solution of the problem. However, additional rigid body 
translations of the rings were necessary in order to insure 
that the axial forces would add up to a zero resultant 
across field B— B* , 

The unbalances were eliminated by displacing points 7, 8, 
and 9 in the x— direction. The required displacements were 
found by solving a 4— by— 4 matrix which included the eaui— 
librium conditions in the x— direction at the three points, 
and the requirement of a zero resultant force in the x— 
direction in a transverse section across field B— B 1 , These 
displacements, of course, gave rise to unbalanced tangential 
and radial forces and to moments in the plane of ring B— B . 

The unbalances were eliminated by undertaking suitable tan- 
gential and radial displacements and rotations at points 7 
and 8 of ring B— B , and a suitable radial displacement at 
point 9, Because of the symmetry point 9 could not 'undergo 
any rotation or tangential displacement. The magnitudes of 
the displacements were determined by solving a 7— by— 7 matrix. 
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The displacements r enuired for balancing the forces and 
moments in the plane of ring B— B threw "back unbalances into 
the x— dir ect ion . They were again eliminated by displacing 
points 7, 8, and 9 in the x— direction and undertaking a 
suitable amount of rigid body translation of ring B— B in 
the x— direction. After this, it was again found necessary 
to balance the plane of ring B— B as before. The displace- 
ments needed for this last balancing caused insignificant 
unbalances (about l/2 of 1 percent) in the x— direction. 

After all these displacements were undertaken, points 
7, 8, and 9 could be considered to be in equilibrium for 
practical purposes. However, a check calculation showed 
that point 6 was considerably out of balance in the plane 
of ring B— B. Hence point 6 was now moved and the un- 
balances thrown back by these motions upon point 7 were 
balanced by moving points 7, 8, and 9 once more in the 
plane of ring B— B. Again sizable residuals were found to 
exist at point 6. The motions of point 6 required for 
balancing the forces and moments at point 6 caused rela- 
tively large unbalances at point 7, It was found that the 
convergence of the procedure could be accelerated by moving 
points 6 and 7 now simultaneously. After this step all the 
unbalances had values which could be considered as negligible. 
The establishment of a check table, however, showed that the 

^AB — &bB * 

resultant force quantities X^.g , ^BB' ’ an ^ ~ 

as well as the residuals at point 5, were too large. A few 
additional motions sufficed for reducing these quantities to 
permissible values. 


Numerical Calculation of the Equilibrium 
of the Small Cylinder 

The complete operations table is presented as table 2. 
Because of its large size it is shown symbol ically , and the 
symbols are explained in table 3, Symbols P, Q, R, and S 
have two values each, one corresponding to the complete 
cylinder (no cutout) and the other to the cylinder with the 
cutout. It is noted that the effective width of sheet attached 
to stringer 8 is reduced when the sheet in the panel between 
stringers 8 and 9 is cut out. The entries in the table 
were calculated according to the principles discussed in the 
first part of this report. 
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Individual displacements of the points situated along 
ring A— A are assumed in spite of the fact that the ring has 
to he considered as a rigid body. In the establishment of 
the operations table symmetric motions of the entire system 
were considered throughout. Consequently, whenever the 
nth point of ring B— B is moved through a unit distance, say 
in the tangential direction, point n' of ring B— B , point n 
of ring B'~ B T , and point n' of ring B' also are moved 
simultaneously likewise, The effect of these simultaneous 
motions was duly considered when the operations table was 
established. The effect is noticeable in the entrieB 
referring to stringers 3 and 8. The points on ring B— B 
also are affected. The unit displacement is 0.001 inch 
and the unit rotation is 0.001 radian. 

Similarly, whenever a point on ring A— A is displaced, 
the three points symmetrically situated to it are also dis- 
placed. A motion of a point on ring A*— A* , however, has 
no effect upon the forces listed in table 2. On the other 
hand, the displacement of a point on ring A— A and stringers 
3* and 8* influences the entries in rows 2 and 8 in table 

2, as may be seen from the numerical values given in table 

3. 


Table 4 lists the forces and moments caused by rigid 
body displacements. A rigid body displacement is defined 
as a set of displacements of a number of points during 
which the distances between the points do not change. 

The effect of a rigid body displacement can always be calcu- 
lated as the sum of the effect of the individual displace- 
ments that constitute the rigid body displacement. In table 
4 all values are listed for the complete cylinder (no cutout). 
The rows marked contain the forces and moments caused 

by a rigid body rotation of ring A— A , the rows marked u>-g-g 


those caused by a rigid body rotation of ring B— B , and the 
rows marked T|^ those caused by a vertical downward transla- 
tion of ring BB , The magnitudes of these motions are so de- 
fined that the displacement of point A1 is —0.003 inch in 
the x— direction for the rows marked and the displace- 

ment of point B x is —0.001 inch for the rows marked id^. 


In both cases 
of rotation. 


the horizontal diameter of the ring is the axis 
In the downward translation corresponding to 


the rows marked the displacement is 0.001283 inch 


It is noted that a vertical downward displacement n of 
the intersection point of a stringer with a ring defined by 
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an angle <p measured from the vertical direction as shown 
in figure 27 can he considered as the sum of a tangential 
displacement of sintp and a radial displacement of p 

coscp. The sum total , of all the rigid "body displacements 
listed is a group operation which gives vanishing residual 
forces and moments throughout except for the forces acting 
on ring A— A in the x— direction. The latter add up to a 
pure counle of 49,000 inch— pounds about a horizontal trans- 
verse axis . 

In some instances consideration of all the displace- 
ments simultaneously reduces the work of computation. ffor 
instance, in the case of a rigid body displacement of a ring 
In the x— direction it is self-evident that no shear will 
arise in the panels. The only forces caused by such a 
displacement are due to the shortening of the stringers. 

The displacements contained in table 4 constitute the 
solution of the problem of bending of the cylinder not having 
a cutout. Since it is anticipated that the effect of the 
cutout will be restricted to the immediate neighborhood of 
the cutout, in the balancing procedure to follow only a 
portion of the operations table will be used. This por- 
tion is presented as table 5. It contains only the displace- 
ments of points 7, 8, and 9 and four rigid body translations 
of the rings, all related to the cylinder with the cutout. 

The upper— left 10— by— 10 corner of the operations table 
is identical with the corresponding portion of table 2: 


M L 
n p q 
a s 

The eleventh row represents a rigid body translation 
In the x— direction of ring A— A. The eleventh column contains 
the contributions of the individual and group displacements 
to the resultant axial force acting In a complete transverse 
section of field A — B. Similarly, row 12 corresponds" to a 
rigid body x— translat ion of ring B— B (combined with a 
simultaneous symmetric translation of ring B f — B ' ) , and column 
12 contains the resultant x— forces in a transverse section of 
field B— B 1 , K.ow 13 contains the f orces and moments caused by 
a rigid body rotation of ring B— B about its horizontal 
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diameter such that i — —0.001 inch for the intersection 
point of stringer 1 with ring B— B . Simultaneously, of 
course, ring B l — B* is also rotated symmetrically. In 
column 1?- are listed the contributions of the individual 
and group operations to the expression ( 

where M^ B is the moment about the horizontal diameter of 

the cylinder in a transverse section through field A— B , 

^BB corresponding quantity in field B— B* , and r is 

the radius of the cylinder. Finally the last row represents 
a rigid body vertical downward displacement of ring AA , 

(Of course ring A’— A' must also be displaced simultaneously, 
but this displacement will not have any effect upon the 
quantities listed in the operations table.) The magnitude 
of the translation is such that the radial displacement v 
of the intersection point of stringer 1 with ring A— A ig 
+0.001 inch. The last column contains the contributions 
of the individual and group displacements to the vertical 
shear force resultant acting on ring A— A. 

On the assumption that the displacements calculated in 
table 4 for the complete cylinder represent in first approxi 
nation the displacements of the cylinder with the cutout, 
the values are substituted in the operations table of the 
cylinder with the cutout. The only unbalances corre- 
sponding to these displacements occur at points B8 and B9 
in the x— direction, because the displacement pattern of 
table 4 does not contain shearing deformations in field 
BB 1 . The unbalances are 


X B8 = pounds X Bg * 613.32 pounds (25) 

Oorr esponding to these values the sum of the axial forces 
in a section through field B— B* is not zero but 

X BB' = 755 *6? pounds (26) 

and the difference in the moments in the sections through 
fields AB and BB* is (divided by the radius of the cylinder) 

)/r = 744.8 inch— pounds (2?) 

Points B7 , B8 , and B9 are now balanced, and the axial 
force resultant in field B— B* is reduced to zero by suitable 
x— displacem ent s of points B? , B8 , and B9 , and by suitable 
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rigid "body translation* of ring* BB and AA in the x— direction. 
The upper left 3— by— 3 corner of the matrix shown in table 
6 contains the axial forest caused by the axial displace— 
aents of the three point* in question* The fourth row 
gives the axial forces caused by a unit axial rigid body 
x-trenslat ion of ring B— B combined with an axial trans- 
lation of ring A— A through a distance of three units in 
the x— direction. Bing A— A has to be moved in order to in- 
sure that points Bl to 56 be net thrown out of balance be- 
cause of the rigid body translation of ring B— B. The com- 
bined operation thus defined is denoted as | group. It 
may be seen that all the figures in table 6 either are taken 
directly from table (5 t Of are combinations of values listed 
in table 6* 

In order to balanoe the xi-resi&uals obtained, the eet 
of equations represented by the matrix of table 6 and the 
right-hand side members given in equations (25) and (26) 
are solved by the matrix method. The results are 


l 7 .= 0.0253 $ 8 * 0.12816 | 9 - 1.9190 £g roup * 0.07 (28) 

The values are given in 1/1000 inch as usual. 


Substitution of the aboTe displacement values into part 
of table 5 gives the residual forces and moments acting in 
the plane of ring B— B at the location of stringers 7, 8, 
and 9. (The residual* in the x— direction were balanced 


out in the preceding step of th 
residuals are; 

T p7 * —0,480 pound 
T B8 * —21.868 pounds 
H B7 * ~° .07266 pound 
H-gg » —1.5798 pounds 


calculations.) The new 


* 3.352 pounds 
N B7 « —0.03717 inch— pound 

( 29) 

Ngg * —0,5487 inch— pound 


To eliminate the unbalances listed in equations (29) the 
seven equations represented by the matrix of table 7 (which 
is just another part of table 5) together with the right- 
hand members given in equation* (29) are solved by the 
matrix method. The solution let 
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U B7 " °* 079 

V B? 

= -0.850 

w B? = -1,429 


Ujjq = —1.079 

V B8 

= -2,602 

W B8 s ^* 9 ®4 

(30) 


V B9 

= 6.68 




When these displacements and. rotations are undertaken, the 
forces and moments acting in the plane of ring B— B are in 
equilibrium at points B7 , B8 , and B9 . However, the equi- 
librium of the forces acting at these points in the x-direc— 
tion has been disturbed. The unbalances thrown back in the 
x— direction are calculated again from table 5. They are 


= --11,06 pounds X B8 = —10.86 pounds = 41.18? pounds 

(31) 


No unbalaneed axial force results in a transverse section 
of field BB * . The residuals are small as compared to the 
original ones. Nevertheless, they are eliminated by using 
once more the matrix of table 6. The necessary displacements 
are 


^B7 

* -0.01174 

l Bg » 0.1232 

^B8 

= -0.00S684 

1 m 0.00274 

group 


(32) 


After these axial displacements were undertaken, the unbalances 
in the plane of the ring are: 


T b? = 0.222 

* 

B7 

-0.006 

N b? = 0.0188 


T Bg = -1.604 

R B8 = 

-0.062 

N Bg = —0,054 

(33) 


R B9 " 

0.371' 



Use of the matrix of 
ments and rotations: 

table 

7 gives the 

following displace— 

U B7 = 0,010 

t B7 

* -0.0637 

w B? » —0.121 


U B8 ' ■’-°» 0844 

V B8 

* -0.213 

W B8 = 0,144 

(34) 


V B9 

■ 0.55? 
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The axial unbalances caused "by these distortions are negli- 
gibly small. 

Substitution of all the preceding displacement "values 
into the operations table (table 5) reveals that a differ- 
ence now exists between the moments transmitted through 
fields A— B and B— B': 


( K AB — M BB*)/ r = ~ 28 * 78 inch— pounds per inch (35) 

This can be eliminated by rotating ring B— B through an angle 
( in 1/10,000 rad) 


U) SB = —0.005 


(36) 


The unbalances caused at the different points by this rotation 
are found to be negligibly small. 


Points B7 , B8 , and B9 can now be considered as completely 
balanced. Substitution of the displacement values corre- 
sponding to all the individual displacements of point B7 in- 
to that portion *f table 5 which represents the int er 1 ink: age 
between points B6 and B7 shows that point B6 is out of balance. 
The unbalances are 


Xpg = 2.082 pounds 

T = 3.01 pounds R__ = 1.123 pounds N 

DO DO 


(37) 

0.893 inch— pound 


The residual forces and moment listed in equations (37), 
together with the matrix of table 8, constitute a system of 
four linear equations which permits the calculation of the four 
displacements of point B6 necessary to balance out the r « * idual s . 
The displacements are 

£ = 0.002582 

B6 

U B6 = -0.01287 v B6 = 0.20924 w Bg = 0.28929 (38) 

These displacements, while balancing point B6, throw unbalances 
upon points B5 and B7, The former were not recorded: the 
lat t er f ollow; 
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T_ = 4.1868 pounds 
B7 

X B7 = —0.0826 pound 


R^ 7 = —0,9831 pound 


N 


B7 = 0.3 609 inch— poundj 


>( 39 ) 


The three quantities pertaining to the plane of the ring are 
now "balanced out using the matrix of table ?. The displace— 


Bents obtained are: 





u B7 = 0.03607 

t B7 

= -0.23349 

w B7 = —0.11682 


= -0.02859 

V B8 

= 0.01284 

w = 0.07393 

B8 

.(40) 


V B9 

= 0.08032 

/ 



The unbalances caused at B6 by the above— listed displacements 
of B7 are: 


T =3.778 pounds 
B6 

= 0,601 pound. 


R = 0,812 
B6 

M Bg = 0.522 


pound 
inch— pound 


(41) 


These residuals are again eliminated with the aid of table 
8. The following displacements are obtained: 


ujjg = 0,02233 


! B6 = 0.00092 
v B g = —0.15102 


W B6 



The effect of these motions on point 7 is found to be: 


T = 3.33 
B7 


X = -0.3932 
B7 

R = -0.6636 
B7 


N 

B6 


0.402 


>(43) 


Comparison with the values shown in equation (39) in- 
dicates that this process is very slowly convergent, if at 
all. A rapid elimination of the residuals at both points 
6 and 7 can be had only by moving both these points at the 
same time. The motions undertaken are: 
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U B6 

= 0.012 

U B7 

0.019 


V B6 

= 0.11 

V B7 

-0.17 

> (44) 

W B6 

CO 

o 

. 

o 

II 

< 

w 

il 

0.08 

-J 



At this stage of the relaxations all the remaining unbalances 
on points 6, 7, 8, and 9 are considered as negligibly small. 

A check table is set up and is presented as table 9. It 
indicates that point 5 is out of balance. 

The check table also shows unbalances for ^BB* * 

and ( M^ — ^BB t )/ r . ^^ ese residuals probably are due to 

some errors in the numerical calculations. All the residuals 
are reduced to negligibly small quantities by additional 
operations contained in table 9. 

Figures 28, 29, and 30 contain the axial stress, the 
bending stress, and the shear stress distributions, respec- 
tively, in the small cylinder as calculated from all the 
displacements determined in this section. 


NUMERICAL CALCULATION OF THE EQUILIBRIUM 
OF THE LARGE CYLINDER 


It was hoped that application of the procedure just 
shown would result in establishing the equilibrium of the 
large cylinder in a reasonable number of steps. This an- 
ticipation, however, was not fulfilled and the calculations 
became so time consuming that they cannot be recommended 
for routine work, although the results obtained were in 
good agreement with experiment. 

The system of designating the individual points is shown 
in figure 31. As may be seen, the "large’* cylinder is' just 
a non— s impl if ied version of the same cylinder that was calcu- 
lated by the matrix method in the first part of this report. 
Only one— quarter of it need be considered because of the 
symmetry. This quarter contains one rigid and three non- 
rigid half— rings (one of the latter is cut) with altogether 
35 points which have a total of 97 degrees of freedom of 
motion. 
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(Che operations table is represented symbolically in 
table 10. The squares denoted by 1, 2, and 3 are identical 
with the squares that are similarly situated in table 2. 

Table 11 contains the square designated 4, and table 12 
those denoted by 5, 6, 7, and 8. Symbols A to Z in these 
tables have the same meanings as before (see table 3)* and 
the symbols a to u are explained in table 13. 

The principle used in solving this operations table was 
the same as that discussed in connection with the small cyl— 
linder. First rigid body displacements were undertaken with 
all the eight rings in order to find the solution for the 
complete cylinder (no cutout). This involved rotations of 
rings D— D , C— C , B— B , and A— A in the ratios 1:3:5:7, and 
vertical downward translations of the last three in the 
ratios 1:3.02:6.01. Next the displacements obtained were 
substituted in the operations table for the cylinder with 
the cutout and the unbalances were calculated. These were 
then balanced out by solving the matrix of all the x— dis- 
placements considering only points 6 to 9 on rings B— B and 
0— C , and 5 to 8 on ring D— D . The displacements undertaken 
caused unbalances to arise in the planes of the three rings. 
Three matrices were set up to take care of the forces and 
moments in the plane of each ring individually. The un- 
balances in the plane of one of the rings were eliminated 
first by solving the corresponding matrix. The unbalances 
caused by these displacements in the plane of the next ring, 
together with the original unbalances there, were then 
balanced by solving the corresponding matrix, and so on. 
Because of the unexpectedly strong interaction between the 
rings the matrices had to be solved many times before the 
unbalances were reduced in all the rings simultaneously. 

The displacements undertaken in the planes of the rings 
during this balancing procedure three unbalances back in 
the x— direction which necessitated a repetition of the 
entire procedure. 

At the beginning, the unbalances in the x— direction de- 
creased after each complete balancing in the plane of the 
rings but later they began to increase gradually. At the 
same time, the displacements in the neighborhood of the 
cutout increased steadily and tended to attain unexpectedly 
large values. Consequently, the procedure adopted was found 
to be divergent. It is possible, however, that the divergence 
was caused either partially or wholly by a slight error In 
the operations table. 
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It was also observed that point D8. underwent large out- 
ward displacements (in the negative r direction) while in 
the solution of the simplified structure in the first part 
of this report the displacement of the corresponding point 
was small and inward. Because of this the displacement 
pattern of the rings was arbitrarily changed to conform 
better with that found in the case of the simplified 
structure, The procedure of balancing was then continued 
as before and was found to converge, though slowly. It 
might be mentioned that in the tests described in reference 
6 both inward and outward deflections were observed. 

At suitable stages of the procedure again some additional 
points (B5, 05, 04, and D4) had to be displaced in order to 
reduce the? residuals all over the structure to negligibly 
small quantities. The final displacements obtained are 
listed in table 14, and the final residual forces and moments 
in table 15. The stresses calculated from the displacements 
are shown in figures 19, 20, 24, and 26. 

It may he seen from figures 19 and 20 that the axial 
stress distribution is very much the same in the solutions 
corresponding to the simplified cylinder and the large 
cylinder. The straight-line portions of the diagrams are 
practically parallel although not coincident. The reason 
for the shift is that the location of the centroid of the 
cross section of the simplified cylinder Is not the same 
as that of the large, and consequently also the actual, 
cylinder. It can be anticipated, therefore, that the agree- 
ment is better between the experimental curves and the 
theoretical curves calculated for the large cylinder than 
between the experimental curves and the theoretical curves 
caleulated for the simplified cylinder. This is borne out 
by figures 19 and 20. Altogether the agreement between 
theory and experiment is good. 

In the shear curves of figure 24 and the bending moment 
curves of figure 26 the agreement is good between values 
caleulated for the complete portion of the cylinder on the 
basis of the simplified and the large cylinders. Consider- 
able deviations occur in the cut field. This could be ex- 
pected since the simplifying assumptions changed the mechani- 
cal conditions in this region. 



NACA TN No, 1014 


34 


A SIMPLIFIED APPROXIMATE SOLUTION 


The problem can be simplified radically by assuming 
that the rings are infinitely rigid in their planes and by 
establishing the equilibrium of the i— forces only. This 
proposition was worked out on the basis of the operations 
table of table 1. The elements related to the equilibrium 
of the axial forces are contained in the upper— left 9— by— 9 
corner of the operations table. 

The equations corresponding to the assumed rotation 
of the rigid end ring were solved by the matrix method. 

The displacements of points A, B, 0, D, and E were found to 
be 


C A ■ -0.1482 t - -0.4340 


f = -0.3799 

V 


= -0.0060 


! e = 0.5488 


> (45) 


The axial stresses were calculated from these displace- 
ments. They are shown in figures 33 and 33 together with 
the curves of axial stress calculated in the first part of 
this report. The agreement was found to be excellent between 
the present approximate solution and the exact solution of 
the problem of the simplified cylinder. Of course, the 
approximate solution does not give any useful data for the 
calculation of the bending moments in the rings and the shear 
stresses in the sheet covering. 


00N0LUS IONS 


The stress distribution caused by a pure bending moment 
in a cylindrical reinforced monocoque cylinder having a 
symmetric cutout was investigated by several methods of calcu- 
lation. The results were compared with data obtained in ex- 
periments described in reference 6. The main conclusions 
follows 

1. The operations table of the problem as defined in the 
Southwell method can be set up easily if use is made of the 
formulas contained in figures 7 to 10. 
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3. The set of simultaneous linear equations represented 
by the operations table can be solved by the matrix method 
as shown in reference 4 if the number of unknowns is not 
too great. (The 30— 'by-»-30 matrix shown in table 1 ean be 
solved in from 1 to 4 days depending upon the operator and 
the calculating machine.) 


3. The calculated axial stresses are in good agreement 
with the experimental data presented in reference 6. (See 
figs. 19 and 20.) The difference in the location of the 
neutral axes corresponding to test and calculation (the 
latter labeled simplified cylinder) is due to the fact that 
the location of the centroid of the actual structure differs 
from that of the simplified structure. 


4. The step— by— step procedure developed in the second 
part for solving the operations table of the so-called large 
cylinder (table 10) was slowly convergent and is not recom- 
mended in its present form for practical use, The results 
obtained by it for the axial stresses are in good agreement 
with test results. (See figs. 19 and 20.) 

5. Sxperimental values for the shear stress in the 
cylinder shown in figure 1 were available only in the cut 
section. They do not compare favorably with the calculated 
values. However, the experimental values are not considered 
reliable as stated in reference 6. 


Polytechnic Institute 
Brooklyn , N . Y , , 


of Brooklyn, 
duly 1945. 
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APPENDIX 


To show how figures 7 to 10a representing the four panel 
problem may be put to use in setting up the operations table, 
two numerical examples are worked out. Pirst , use is made 
of figure 7 to determine the forces and moments introduced 
at the constraints when point B of the simplified structure 
is displaced axially through a positive unit distance of 
0.001 inch. The following data are needed in the f 1 aleU'r» 
lations: 


E . = 

str 

10.3 X 10 s 

pa i 


& 

strpp 

■^etrpH 

0.1877 

in.® 

a II 1 IV 

= 12.86 in 

• 


a X,III 

«* 9 . 64 in . 



*I,IX 

* L III, IV 

3 .927 

in . 

G II " 

G 1II = & IV 

= 0.385 

E 4. 
str 


G-j = 0 

tjj = tjjj =• t jy .= 0.312 in. 
(Ot)| = 47.64 lb. 




( la) 


3 .97 'X' 10 6 -pui 


J 


The values of the coefficients a, were calculated from 
the simplified formulas suggested in the conclusion of refer- 
ence 5. They were checked by the values taken from figures 
86 to 93 of reference 5 for the smallest and largest values 
of | and V, respectively. All calculations were carried 
out by slide rule. The results are: 
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(a t ) I, 

XI 

“ ^ ^ 1 1 1 , IT 

= -0.499 




(a r ) I, 

II 

= ^ a r^IH,IT 

= 0.039 


> 

(lb) 

(*n^I f 

II 

= < a n ^ 1 1 1 , IT 

* -0.0033 





For A the forces and moments acting upon the constraints 
then become, according to figures 7 and 10a: 


X A = (47.64)( 12.86)/( 4)(3 .987) = 38.9 lh 


T a = -(0.499 /2) 47.64 
S A = -(a 039/2) 47.64 


= -11.9 lb 
= -0.934 lb 


>( 2a ) 


ff A = _(o. 0033 /2)( 47.64 )( 3.9 27) = -0.305 in.-lbj 


At E the motion causes: 

3 


X B = -10.3 X 10 (0.1877)[(l/l2.86)+(l/9.64)] 


+(47.64)[2(l2.86) + 9.64'Tl/( 4)(3.927)] =-433 lb 


T B = -( 0.499/2)47.64 


= -11.9 lb 

R b = (a039/2) 47.64 = 0.934 lb 

N b = -(0.0033/2)(47.64)( 3.927) = -0.305 in. -lb 

At point 0 there results: 

X Q = ( 47. 64)( 12.86 + 9 .64) /( 4)( 3 .9 37> - 68.1 lb 

*0 - R c - N 0 - 0 


( 2b) 


( 2c) 


0 
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For 3? and G the values are: 

Xp a [10.3 X 10 3 ( 0.1877)/( 9 .64) ] - ( 47.64)( 9 . 64) [l/( 4)( 3 .927 ) ] 
= 146.8 lb. 

X = ( 47.64)(9 .64)/( 4)( 3 .927) = 29.2 lb 
G 

T^ = T & = 11.9 lb 

= 0.934 lb (2d) 

Hj, = R & a 0.305 in .—lb 

To illustrate further the use of figures 7 to 10a, the 
effect of a tangent ia.l displacement of point G through a 
unit distance of 0.001 in. is investigated with the aid of 
f igur e 8 . 

In the course of finding the ring influence coefficients 
needed in the calculations, the ratios 

(EI) r /L, (EI) r /L 2 , (EI) r /L 3 

have to be determined. The moment of inertia of the ring 

plus its effective sheet is found to be 8.05 x 10 s in. 4 ; 

Lj is 3.927 in. for arc EG; and Ljjj jy = 7.854 in. for 

ring segment GH. 

Further, convenient values of the parameters Y and £ 
have to be assumed before use can be made of figures 14 to 
31 for the movable end influence coefficients and figures 
50 to 67 for the fixed end influence coefficients, or tables 
III and IV, all of reference 5, The values Y = 10,000 and 
I = 0.25 were found to be the closest choice for the given 
ring elements. For FG, which subtends an angle of 22.5°, 
the graphs were used; for GH the necessary data were taken 
from the tables. The final results are: 
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Jor arc ^G: 


nnj^ s 1.539 

nny 

= 

0.2853 

A 

£ry ** —1 .378 

nry 

= 

-0.792 


n ^M “ 5*38 

nt j> 

= 

in 

« 

in 

> 

= 1 * 777 

”» 

= 

1.461 


tr^ * —7 .95 

try 

s 

-8.10 



1 tjl 

= 

41.4 

J 

For arc GH: 

nn M = 0.874 

nn p 

- 

0.3562 

"\ 

nrjj - —0.423 

. — -*» 

nr y 

s 

-0.273 


n ^M ~ 

^ty 

3= 

0.871 


rr„ =« 0 . 28 1 

rr F 

s 

0.2445 

► 

tr M » -0,6275 

t7y 

rc 

-0.641 


tty ” 1.538 

ttjl 

ss 

1.530 

J 


> (3a) 


> (3b) 


These Talues correspond, to a unit displacement of 0.001 
in. or a unit rotation of 0.001 rad. 

When point G is moved, Gj and Gjy set equal 

to zero since panels I and IV according to the notation of 
figure 8 of this report are cut out. The parameters 

, a r , a n for arc FG are identical with those listed in 

equation (lb). Those for arc GH were derived in analo- 
gous manner and are: 



Pi rt 
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at 1 1 1 , 17 


= -0.490 


a riH,IV = 0,079 

a = - 0.0065 

n III, IV 


f (Sc) 


J 


The remaining geometric and mechanical properties occurring 
in the calculations are identical with those given by- 
equation (la). 


The forces and moments arising at the constraints when 
oint & is moved tangentially threugh the unit distance can 
ow be written down. 


X B = (0.499/2) 47.64 
T b = (0.499) S (47.64)(3.927)/9.64 
R b « (0.039)(0i499)(47.64)(S.927)/9.64 
N b * (QL0033)(0i499)(47.64)(3.927) 2 /9,64 


* 11,9 lb 
= 4.85 lb 


- 

> (4a) 


«= 0,378 lb 

=* 0 1 126 in.-lb^ 


Xq = -( 47.64)( 0.499 — 0.49 0 )/2 ’ = -0.1 lb ^ 

Tq = 47. 64[( 0,499 ) S ( 3 . 927 )+( 0 . 49 0 ) ®( 7 . 854)]/9.64 = 14.16 lb 

R 0 - 47. 64[( 0.490)(0.079 )(7.854)-( 0.499)(0.039 )(3, 927)3/9. 64 >( 4b) 

- 1.126 lb 

H 0 = 47. 64[0. 490)(0.0065)(7. 854) a +( 0.499)(0.0033)(3. 927)]/9. 64 


1.099 in .-lb 
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Xjj = -( 0.49 0/2) 47.64. 

T d = ( 0.490} S ( 47.64)( 7.854)/9 .64 

R d = -(a079)( 0.490)( 47.64)(7.854)/9.64 

N d = (a0066)( 0.490)( 47.64)( 7.854) */9 .64 


= - 11.8 11 

= 9.31 11 


= -1.504 11 


= 0.973 in. -11 

-j 


x r = x b “ 11 * 9 1V 

= 41.4 - E = 36.55 11 
Ry = 8.10 - R B = 7.722 11 
Ny = 5.54 — Ng = 5.414 in.— 11 


> 

- 


Xg_ *= Xq = —0,1 11 *> 

= -41 -1.533 - T 0 - - 56.698 11 
R g = 7.95 -0.6275 - R c « 6.1965 11 
N 0 = -5.38 -0.81 - J5 = -7.289 in. -11 

w C 

X H " X 0 * “11-8 1* 

T = 1.530 - T = -7.78 11 

8 B } 
R H = -0.641 -R d = 0.863 11 

N h = 0.871 -Njj = -0.102 in .—11 


^ ( 4c) 


( 4d) 


( 4e) 


(4f) 
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TABLE I. OPERATIONS TABLE FOR THE SIMPLIFIED CYLINDER. 
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Table H. Complete Operations Table for Small Cylinder 


Stringer No. 


_ 1 

2 

3 

4 

Ring AA 5 
6 

7 

8 
9 


Ring AA 

123456789 
A C 
BED 
FED 
FED 
FED 
FED 

FED 
F E G 


Ring BB 

123456789 
T V 
U X W 
Y X W 
Y X W 
Y X W 
Y X W 

Y X W 
Y X Z 
V X T 


1 

2 

3 

4 

Ring BB 5 
6 

7 

8 
9 


T V» 

V»X W» 

Y»X W« 

Y»X W« 

Y‘X W« 

Y»X W» 

Y’X W» 
Y*X Z‘ 
VjT 


H K 
JHL 
N M L 
N M L 
N M L 
NML 
N M L 
N P Q 
R S 



Table TTT. Sections of Operations Table for Small ^Cylinder 



I 

T 


R 

N 


X 

T R 

pi 


I 

T 

pi 

N 


Section f\ 






Section F 




Section 

1 



-344-5 



-1.872 


* 

18.91 

11.55 0.936 

0.29 

1 

18. 91 

11.55 

-0.936 

0.29 

V 

-1.872 



-5.39 


u 

-11.55 

53.45 H.97 

8.09 

u 

-11.55 

53.45^ 

-11. 97 

8.09 







V 

0.936 

-11.97 -2.28 

-1.30 

V 

-0.936 

11.97 

-2.28 

1.30 



Section 

B 



yt 

-O.29 

8.09 1.30 

O .63 

w 

-0.29 

8.09 

-1.30 

0.63 

Y 

37.82 



1.872 


* 


Section 0 



Section U = P(no cutout) 

u 

-23.1 



23.94 











V 

1.872 



-4.56 



37.82 

1.872 


\ 

-957.92 0 

1.872 

0 

It 

-o;58 



2.60 


u 

23.1 

-23.94 


u 

0 

-139.9 

0 

-16.63 


r 




V 

1.872 

-4.56 


V 

1.872 

0 

-5.39 

0 



Section 

C 

* 



0.58 

—2. 60 


V) 

0 

-16.63 

0 

-3.86 

Y 

18=91 

-11. W 


0.936 

-0^29 


Section H ■ ,s(no cutout) 



Section 

N 


V 

0.936 

11.97 


—2.28 

1.30 








-0.936 

-0.29 






'n 

-957.92 1.872 


\ 

18.91 

-11.55 



Section 

C 1 



V 

1.872 

-5.39 


u 

11.55 

53.45 

11.97 

8.09 





dr 






V 

-0.936 

-11.97 

-2.28 

“1.30 

1 

18. 91 

+11.55 


0.936 

0.29 



Section J 


it 

0.29 

8.09 

1.30 

0.63 

V 

0.936 

-n.97 


-2.28 

-1.30 








0 cutoutv) 






* 

37-82 

-1.872 



Section Q (n 




Section 

D 



u 

23.I 

23.94 

f 












V 

-1.872 

-4.56 


% 

37*82 


-1.872 



18.91 

-11.55 


0.936 

-0.29 

n 

0.58 

2.60 


..u 

—23.1 


-23.94 


u 

11.55 

53.45 


-11.97 

8.09 





V 

-1.872 


-4.56 


V 

0.936 

11.97 


—2.28 

1.30 



Section K 


n 

-0.58 


-2.60 


11 

0.29 

8.09 


-1.30 

0.63 

















18.91 

+11.55 -0.936 

+0.29 


Section S (cutout) 




Section 

E 


* 

V 

-0,936 

-11.97 +2.28 

-1.30 




1.872 












-344.5 




-344*5 

0 


-1.872 

0 


Section R (no cutout) 

▼ 

1.872 


-5.39 


u 

0 

-139.9 


0 

-16.63 



-11.55 -0.936 







V 

-1.872 

0 


-5.39 

0 


18.91 

-0,29 




♦ 


■W 

0 

-16.63 


0 

-3.865 

▼ 

-0.936 

-11.97 “2.28 

-1.30 
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Table III. Sections of Operations Table for Small Cylinder (cont'd.) 



I T 

R N 


X 

T 

R 

N 


X 

T R N 


Section 

R (cutout) 



Section 

V 1 




Section U* 

'S 

18.91 -11.55 

-0.936 -0.29 

* 

13.91 

-11 « 
— » ✓/ 

~Q Q -aA 

yJm /yv 

-0.29 

s» 

5 

ATI AA 

j(.o2 

1.8?2 

V 

-0.936 -11.97 

-2.28 -1.30 

V 

0.936 

-0.571 

-0.0462 -0.014 

u 

23.1 

1^142 









V 

-1.872 

-0.0925 


Section 

P (cutout) 



Section 

w 


It 

0.58 

0.028 

> 

-880.92 0 

1.872 0 


18.91 

11.55 

-0.936 

0.29 



Section V 

u 

0 -139.9 

0 -16.63 

u 

11.55 

7.05 

-0.571 

0.'177 




V 

1.872 0 

-5.39 0 

V 

0.936 

0.571 

-0.0462 0.014 

* 

18.91 

-11.55 0.936 -0.29 

w 

0 -16.63 

0 -3.86 

w 

0.29 

0.177 

-0.014 

0.0044 

V 

-0.936 

0.571 -0.0462 0.014 


Section 

Q (cutout) 



Section 

X 




Section V£ 

* 

37.82 

-1.872 


268.9 

0 

-1.872 

0 

* 

18.91 

11.55 0.936 0.29 

u 

—23.1 

-23.94 

u 

0 

14.1 

0 

0.354 

V 

-0.936 

-0.571 -0.0462-0.014 

T 

-1.872 

-4*56 

V 

1.872- 

0 

0.0925 

0 




w 

-0.58 , 

-2.59 

It 

0 

0.354 

0 

0.0088 



Section W* 


Section T 



Section 

I 



18.91 

-11.55 0.936 -0.29 









u 

-11.55 

7.05 -O.57I 0.177 

* 

268.9 

-1.872 


18.91 

-11.55 

-0.936 

-0.29 

V 

-0.936 

0.571 -0.0462 0.014 

V 

1.872 

0.0925 

u 

-11.55 

7.05 

0.571 

0.177 

It 

-0.29 

0.177 — 0.014 0.0044 




V 

0.936 

-0.571 

-0.0462 -0.014 



Section Y' 


Section U 

It 

-0.29 

0.177 

0.014 

0.0044 













18.91 

11.55 0.936 0.29 


37.82 

-1.872 



Section 

z 


u 

11.55 

7.05 0.571 0.177 

u 

-23.1 

1.142 






V 

-0.936 

-O.57i-O.O462 -0.014 

V 

1.872 

-0.0925 

* 

37.82 


-1.872 


tl 

0.29 

0.177 0.014 0.0044 

w 

-0.58 

0.028 

u 

-23.1 


1.142 




Section Z* 




V 

1.872 


-0.0925 






Section 7 

It 

-0.58 


0.028 



37.82 

1.872 









u 

-23.1 

-1.142 

? 

18.91 11.55 

-0.936 0.29 






V 

-1.872 

-0.0925 

V 

0.936 0.571 

-0.0462 0.014 






n 

-0.58 

-0.028 
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Table 17. Solution of Small Cylinder* Without Cutout 




Ring AA 


Rigid Body 
Motions 

A 1 

A 2 

♦ 

A 3 

A 4 

A 5 

A 6 

A 7 

A 8 

A 9 

^ AA = 3 

928.8 

858.3 

656.7 

355.5 

X Forces 
0 

-355*5 

-656.7 

-858.3 

V 

-928.8 

U»bb = 1 

-303.8 

-280.7 

-214.8 

-116.3 

0 

116.3 

214.8 

280.7 

303.8 

*1 AA “ 1,2a 3 

- 11.52 

-10.65 

-8.15 

*^ 4.40 

0 

4.40 

8.15 

10.65 

li .52 

Sura 

613.48 

566.95 

433.75 

234.8 

0 

-234.8 

-433.75 

-566.95 

-613.48 

W AA c 3 

0 , 

10 . n 

18.75 

24.51 

T Forces 
26.52 

24.51 

18.75 

10.11 

0 

w BB “ i 

0 

3.37 

6.25 

8.17 

8.84 

8.17 

6.25 

3.37 

0 

*lAA- 1-283 

0 

- 13.54 

- 25.02 

-32.67 

-35,36 

-32.67 

- 25.02 

-13.54 

0 

Sura 

0 

-0.06 

—0 .02 

0.01 

0.00 

0.01 

-0.02 

- 0.06 

0 

w AA " 3 

0.43 

0.40 

0.30 

0.16 

R Forces 

* 

0 

-0-. 16 

-0.30 

-0.40 

-0.43 

W BB * 1 

0.14 

0.13 

0.10 

0.05 

0 

-0.05 

-0.10 

-0.13 

-0.14 

1\AA - 1-283 

-0.57 

-0.53 

- 0.41 

-0.22 

0 • 

0.22 

0.41 

0.53 

0.57 

Sum 

0.00 

0.00 

-0.01 

-0.01 

0 

0.01 

0.01 

0.00 

0.00 

cn 

« 

3 

0 

0.26 

0.47 

0.62 

Moments 

0.67 

0.62 

0.47 

0.26 

0 

i — i 

H 

1 ffl 

m 

3 

0 

0.08 

0.16 

0.21 

0.22 

0.21 

0.16 

0.08 

0 

^ AA 3 1,283 

0 

-0.34 

•^ 3.63 

-0.82 

-0.89 

—0.82 

-O .63 

-0.34 

0 

Sura 

0 

0.00 

0.00 

0.01 

0.00 

0.01 

0.00 

0.00 

0 
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Rigid Body 
Motions 


W AA- 3 
u>bb * 1 
AA = 

Sum 


^AA* 3 
-BB-l 


Table 17. Solution of Snail Cylinder Without Cutout (cant'd 

Ring HB 

d * B J 'E h B ? BO BY b a 
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Table V. Portion of Operations Table for Small Cylinder with Cutout 
Force? in lb. and Moments in in. lb. 


Motions 
In 0.001 
of in. 

or rad. ! 

T 

B7 


V 

V 

T B8 R B8 


V 

X AB I BB I 

(%B 

-“bb^ 

?B7 

-957.92 

0 

1.872 

0 

18. 91 

11.55 -0.936 

0.29 0 

0 

613.4 -1226.8 

-1301.2 12.70 

U37 

0 

-139.9 

0 

-16.63 

-11-55 

53.45 -11.97 

8.09 0 

0 

0 0 

0 38.98 

++ 1 

V B7 

1.672 

0 

-5.39 

0 

-0.936 

11.97 -2.28 

1.30 0 

0 

0 0 

0 -0.628 

*B7 

0 

-16.63 

0 

-3.86 

-0.29 

8.09 ^-1.30 

0.63 0 

0 

0 0 

0 O.978 

)?B8 

18. 91 

-11.55 

-0.936 -0.29 -880.92 

0 1.872 

0 37.82 

-1.872 

613.4 -1072 

-1557-8 16.00 

u B8 

11.55 

53.45 

11.97 

8,09 

0 

-139.9 0 

-16. 63 -22.1 

^3.94 

0 0 

0 21.10 

v B8 

-0.936 

-11.97 

-2.28 

-1.30 

1.872 

0 -5.39 

0 -1.872 

-4.56 

0 0 

0 -0.821 

*B8 

0.29 

8.09 

1.30 

O.63 

0 

-16.63 0 

-3.86 -0.58 

-2.59 

0 0 

0 0.536 

/B9 

0 

0 

0 

0 

18.91 

-11.55 -0.936 -0.29 -344.5 

1.872 

306.7 0 

-306.7 8.98 

V B9 

0 

0 

0 

0 

-0.936 

-11.97 -2.28 

-1.30 1.872 

-5.39 

0 0 

0 -0.444 

iZ 

306.7 

0 

Q 

0 

306.7 

0 0 

0 306.7 

0 

-4907.2 0 

0 0 

Itsb 

-920.1 

0 

0 

0 -843.1 

0 0 

0 -306.7 

0 

4907.2 -9047 

755.7 0 

^EB 

-652.65 -6.25 

0.1008-0.1569 

-781.58 

-3.38 0.1316 -0.0849-309.58 0.1425 

0 755.67 -6616 71.86 

*\AA 

6.35 

19.50 

-0.314 O.49 

8.30 

10.55 -0.41 

0.264 8.98 

-0.445 

0 0 

0 -224.04 
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Table VI. X Matrix 
%7 X B8 ^9 


V 

1b 

19 

1 Group 


-957.92 

18.91 

0 

0 


18.91 

•880.92 

18.91 

77.0 


0 

37.82 

-344.5 

613.4 


^BB* 

-1226.8 

-1072 

0 

-9047 


Table VII. Ring Matrix No. 1 



t B7 

Rb? 

N B7 

T B8 

®B8 

N B9 

Rb 9 

U7 

-139.9 

0 

— 16.63 

53-45 

-11.97 

8.09 

0 . 

V 7 

0 

-5.39 

0 

11.97 

-2.28 

1.30 

0 

w 7 

-16.63 

0 

-3.86 

8.09 

-1.30 

0.63 

0 

u 8 

53.45 

11.97 

8.09 

-139.9 

0 

• -16.63 

-23.94 

v 8 

—11.97 

-2.28 

-1.30 

0 

-5.39 

0 

-4.56 

"8 

8.09 

1.30 

0.63 

-16.63 

0 

-3.86 

-2.59 

V9 

0 

0 

0 

-11.97 

-2.28 

-1.30 

-5.39 



Table 

VIII. Ring Matrix No. 

_z 

X B6 

T B6 

R B6 

%6 

y 7 

-957.92 

0 

1.872 

0 

U 7 

0 

-139.9 

0 

-16. 63 

v 7 

1.872 

0 

-5.39 

0 

YI 7 

0 

-16.63 

0 

-3.83 



f 


Table VIII* Check Table 
Unbalances in lbs. and in. lb. 

< 


Notions 

I B5 

.. T B5 

% 

% 

*86 

T B6 

*B6 

**86 

Besults for 
Complete Cyl 




M 





J 6 0.003505 

0.06628 ‘ 

-O.O4O5 

-0.0033 

-0.0010 

-3.3575 

0 

0.0066 

0 * 

0.021463 

0.2479 

1.1472 

0.2569 

0.1736 

0 

-3.0027 

0 

-0.3569 

y 6 0.47026 

“0*4402 

-5.6290 

# -1.0722 

-0.6113 

0.8803 

0 

-2.5347 

0 

*6 0-35857 

0.1040 

2.9008 

0.4661 

0.2259 

0 

-5.9630 

0 

-1.3841 

>7-0.00921 





-0.17416 

0.10638 

0.0086 

0.002? 

Uy O.I4407 





1.6640 

7.7005 

1.7245 

1.1655 

v r l. 31717 





. 1.2329 

15.7665 

3.0031 

1,7123 

fly— 1 • 58683 
0.11848 





-0.4602 

-12.8375 

-2.0629 

-0.9997 

ug-1.19199 

; vg-2. 80219 
; fl 8 2.12196 









> 9 2.0422 
Ju 0.21622 

66*9261 




66.9281 




Jeb 0.07274 

-66*9261 




-66.9281 




n me 

n 

\j 

n ai 10 

r\ 

r\ nr? 1 1 
U f Wi-L 

i nLL 
JL. fUU 

A AIAA 

A nAAA 
— U*LAJU J 

r\ nm a 
\J*UWXVJ 

v 9 7.31332 









Sum 

-0.02202 

-1.5773 

-0.3525 

-0.2U7 

1.5513 

1.8111 

0,1449 

0.1408 


01 

H 
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Table TUI. Check Table (contM.) 
Unbalances in lbs. and in. lb. 


Motions 

%7 

t B7 


%7 

%8 

t B8 

®B8 

% 

Results for 
Complete Cyl 
^6 0.003505 

• 

0.06628 

0.0405 

-0.0033 

0.0010 

71.17 




0.021463 

-0.2479 

1.1472 

-0.2569 

0.1736 





v 6 0.47026 

-0.4402 

5.6290 

-1.0722 

0.6113 





"6 0.35857 

-0.1040 

2.9008 

-0.4661 

0.2259 





7 7 -0.00921 

8.8224 

0 

-0.01724 

0 

-0.17416 

-0.10638 

0.0086 

-0.0027 

u^ 0.14407 

0 

-20.1554 

0 

-2.3959 

-I.664O 

7.7005 

-1.7245 

1.1655 

V7 -1.31717 

-2.4657 

0 

7-0995 

.0 

1.2329 

-15.7665 

3-0031 

-1.7123 

Vj -1.58683 

0 

26.3890 

0 

6.1252 

0.4602 

-12.8375 

2.0629 

-0.9997 

0.11848 

2.2405 

-1.3684 

-0.1109 

-0.0344 

-104.3714 

0 

0.2218 

0 

us -1.19199 

-13.7675 

-63.7118 

-14.2681 

-9.6432 

0 

16.67594 

0 

19.8228 

v 8 -2.80219 

2.6228 

33.5422 

6.3890 

3.6428 

-5.2457 

0 

15.1038 

0 

n 3 2.12196 

0.a537 

17.1666 

2.7585 

1:3368 

0 

-35-2882 

0 

•*•8.1908 

7 9 2.0422 
0.21822 





38.6180 

-23-5874 

-1.9115 

-0.5922 

66.9281 




66.9281 




y bb °-°™ 

-66.9281 




-61.3271 




WgB — 0.005 

3.2633 

0.0313 

-0.0005 

0.0008 

3.9079 

0.0169 

-0.0007 

0.0004 

v 9 7.31332 





-6.8453 

-87-5404 

-16.6744 

-9.5073 

hJUill 

A £nr 1 
v.uw 

1 Aim 

AtUJU.U 

n neir >4 

n n/00 

0 Aflo; 

_n A/aA 

— w*vt/v 

0.0891 

_0,0163 
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Table VIII. Check Table (cant’d.) 
Unbalances in lbs. and In. lb. 


Motions 

*89 ........ 

®B9 

X AB _ 

X AB' 

W 

y aa 

Results for 





r 


Complete Cyl. 

613.32 



755*67 

744.8 


) 6 0.003505 



2.150 

-4.300 

-4.561 

0.0241 

0.021463 

i 





1.0903 

V£, 0.47026 






-0.1599 

i* 6 0.358566 






0.4590 , 

? 7 -0.00921 



-5.649 

11.299 

11.984 

-0.1170 

uy 0.144073 






5 i 6158 

v 7 -1.31717 






0.8272 

-1.58683 

- 





-1.5519 

0.11048 

4.4809 

-0.2218 

72.673 

-127.01 

-184.56 

1.9667 

Ug -1.191988 

27.5350 

28.5362 




-25. 1510 

v 8 -2.802185 

5.2457 

12.7780 




2. 3006 

w 8 2.12196 

-1.2307 

-5.4959 




1.1374 

7 9 2.0422 

-703.538 

3.823 

626.343 


-626.343 

18. 3390 

0.21322 

66.9281 


-1070.849 




?bb 0.07274 

-22.3094 


356.950 

-658.08 

54.970 


4» b ^.°°5 

1.5479 

-0.0006 


-3.778 

, 33.308 


t 9 7.31332 

13.6905 

— 39.4188 




-3.2471 

Sum 

5.6706” 

0.0001 

-18.382 

-26.199 

29-598 

1.5332 


01 
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Table U. Final Relaxation Table 
Unbalances in lb. and in, lb. 


Motions 

%4 

T B4 

%4 

%4 

*B5 

t B5 

of Table VU 1 





*0.02202 

-1.5773 

Tj * -0.1 

0.0936 

1.1970 

0.228 

0.1300 

-0.1872 

0 

A All 

n*j — — v«vu 

_ A AO*>A 

A L 1 m 
-v t U4 f c 

. lAl 

— .JA ) t * 

« nrrt 1 

ft 

V 

*1 nftft • 

« o.oe j 





-0.0749 

-0.9576 

U7 a 0.01 

W*BB“ 0.004 

1.4128 

0.0327 

0.0002 

0.0008 

0 

0.0354 

Final Residuals 

1.4832 

0.5825 

0.1242 

0.0804 

-0.2037 

-1-2399 

Motions 


*B6 

®B6 

% 

*B7 

t B7 

of Table VIH 

1.5513 

1.8111 

0.3449 

0.1408 

0.6054 

1.6310 

Vj m -0.1 

0.0936 

-1.197 

0.228 

-0.13 



Hr j, _n na 

n.rtfy 

_n.A/70 

v ' '"t f 

n.in/. 

'r' Ti»w»r 

-rLfrcru 



v6 * o.oe 

0.1498 

' 0 

-0.4312 

0 -0;0749 

0.9576 

U7 - 0.01 

0.1155 

0.5345 

0.1197 

0.0809 

0 

-1.399 

WBB“ 0.004 

-1.4128 

-O.0327 

0.0002 

-0.0008 -2.6106 

-O.0250 

Final Residuals 

0.5206 

0.4686 

0.1656 

0.0405 -£.0801 

"I.I446 

Motions 

% 

T B8 

% 

*B8 

*B9 


of Table VUI 

2.6894 

-0.6496 

0.0891 

-0.0163 

5*6700 

0.0001 

Tt7 “ 0.01 

1 

-0.1155 

0.5345 

-0.1197 

0.0809 



^Qroup * *^*002 

-0.154 




-1.226 


* 0.004 

-3.1263 

-0.0135 

0.0006 

-0.0003 

-1.2383 

0.0005 

m — n n t J...1 

PlUttX flPBJLUUftl P 

n frrV! 1 
*V. JWIf 

rt 

•Vf4«V 

n r\n 
-v«uj 

A 

9 oncn 
J»***St 

ft /vn4 

V«VUVAI 
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Table X. Complete Operations Table for Large Cylinder 


Rings 

AA 

BB 

CC 

DD 

AA 

1 

2 



BB 

3 

4 

2 


CG 


3 

5 

6 

DD 



7 

8 


Table XI. Section 4 of Operations Table for Large Cylinder 


Stringer No. 
1 
2 

3 

4 

5 

6 

7 

8 
9 


Ring BB 


1 

a 

b 


2 

c 

e 

f 


d 

e 

f 


Ring BB 
4 5 6 


d 

e 

f 


d 

e 

f 


d 

e 

f 


d 

e 

f 


8 


d 

e 

c. 


8 

a 
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Table for Largo Cylinder 
kRing DD 

123456789 
T 7 

U X W 

Y X W 

Y X W 

Y X W 

Y X W 

Y X W 

Y r a 

t u 

H K 

J M L 

N M L 

N M L 

N U L 

N M L 

N M L 

N m n 


P 9 


^ < (ivtf a «i d •<: a <3 d vrf 
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Table HII. Elements of Soctions of Operations Table for Large Cylinder 


I • 

T ' R 

N 


X 

T 

R 

N 


Element a 



- 

Element 

■ ci . 


-689.0' 

0 



37.82 

0 

0 

0 

0 

-5.49 



0 • 

-11.40 

-2.24 

-1.28 


Element b 



. 

Element 

- g 


75.63 

0 



75.63 


0 


0 

22.8 



0 


22.8 


0 

—4*84 



0 


—4*48 


0 

2.56 



0 


2.56 



Element c 



Element k 

(cutout) 


37.82 

0 0 

0 


18.91 

-11.55 

-0.936 

-0.29 

0 

11.40 -2.24 

1.28 


-0.936 

-11.97 

t 

-2.284 

-1.296 


Element d. 



Element j 

(cutout) 


37*82 

0 0 

0 


37.82 


-1.872 


0 

46.4 -11.40 

7-91 


-23.1 


-23.94 


0 

11.40 -2.24 

1.28 


-1.872 


-4-567 


© 

7.91 -1-28 

0.625 


-0.58 


-2.592 



Element e 



Element h 

(cutout) 


-689.0 

0 0 

0 


-631.6 

11.55 

0.936 

0.29 

0 

-154.0 0 

-16.99 


11.55 

-147 

0.571 

-16.81 

0 

0 -5.49 

0 

. 

0.936 

0.571 

-5.54 

0.014 

0 

-16.99 0 

-3.874 


0.29 

-16.81 

0.014 

-3.88 


Element f 



For no cutout 


37.82 

0 0 

0 



h = e 



0 

46.4 11.40 

7.91 



j « g 



0 

-11.40 -2.24 

-1.28 



k - 



0 

7.91 1.28 

0.625 



je - a 




Element i Element r (cutout) 



-344.5 1.872 


¥ 

249.3 -11.55 

0.936 

-0.29 

T 

1.872 -5.39 


u 

11.55 7.05 

-0.571 

0.177 




V 

-0.936 -0.571 

0.0462 

-0.014 


Element m (cutout) 


7 ) 

0.29 0.177 

-0.014 

0.0044 

> 

-823.6 -11.55 0.936 

-0.29 


For no cutout 


u 

-11.55 -70.0 - 13.05 

-8.32 





T 

0.936 13.05 -2.698 

1.972 


r ■= X 

r* «= X 


■w 

-0.29 -8.32 1.972 

-1.932 


s = Z 
t = v x 

s* = Z* 

t» - v£ 



For Cutout 



u => T 




s « t = u = 0 



Element r’ 

(cutout) 



S* = t> » 0 
r? * T) * O — 0 


■> 

249-3 11.55 

-0.936 

0.29 


H v 


u 

-11.55 7.05 

-0.871 

0.177 




V 

0.936 -0.571 

0.0462 -0.014 




71 . 

-0.29 0.177 

-0.014 

0.0044 



w W'vff € 


Table XIV . Final Displacements for Large Cylinder 


Axial Displacements 
in ,001 in. 


01 

03 


ring 

5 

6 
7 

*9 


u 4 

v 4 

"4 


n? 


5 
5 

T 6 

*6 

5 

u 8 

v 8 

"8 


Ring A 

Ring B 

Ring C 

Ring D 

7 

5 

3 ! 

. 1 

0.45828 

0.36236 

0.19605 

0,05024 

0.02706 


-0.07831 

-0.07749 

-0.03364 


0.02165 

0.10335 

0.08315 


0.22974 ! 

0.47192 

0,18628 


I.40256 

3.31413 


Displacements in Planes 

of Rings in 

.001 in, or .001 

7.6981 

3.849 

1.2830 



0 

-0.05020 

0.01 


0.05 

0.36260 

-0.05 


0.04 

0.66364 

-0.12 


0.03615 

0.37441 

0.37385 


0.11458 

0.55414 

2.39213 


0.03Q33 

-1.01502 

0.59187 


0.03642 

-0.22174 

0.51767 


-0.45799 

-2.20055 

-5.37900 


-O.47259 

0.99654 

-4-4485 


-0.36193 

-0.32158 

-3.23801 


-0.63495 

-0.58833 

-4.54688 


0.76479 

-2.15896 

8.41864 


-03405 

-2.1623 

5.32914 


1.73885 

-4.77541 

53.20725 


-0-47318 

3.30103 

17.21916 


| 

-1.58075 

13.36862 



0 
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Table IV. Residuals for Largo Cylinder 
Axial Forces in lb. 


B6 

B7 

B8 

B9 

c6 

C7 


08 

C9 

1.166 

-11.687 

-8.616 

13.493 

6.3U 

-5.775 


-0.514 

' 11.015 

D5 

D6 

D7 

D8 

X AB 



^CD 

%>« 

-5.28? 

-10.412 

-3.088 

13.169 

21.47 

9.36 


— 8.76 

0.01 



Ring 

B - Forces in lb. 

Moments in 

in 

. lb. 


T 4 

*4 

*4 

T 5 

% 

• "5 


% 

*6 

1.4014 

0.1172 

0.0403 

1.8386 

-0.0096 

0.0381 


0.5852 

-0.1127 



Ring 

C - Forces in lb. 

Moments in 

in, -lb. 


1.3795 

-0.1182 

-0.0509 

0.7087 

0.0572 

-0.0440 


0.1126 

-0.0384 



Ring 

D - Forces in lb. 

Moments in 

in, 

.-lb. 


0.1882 

0.0241 

-0.0143 

-1.0753 

0.1361 

0.0219 


-0.0938 

-O.OI 96 



Ring 

B - Forces in lb. 

Moments in 

in. 

,-lb. 


«6 

®7 

R? 

H 

. t 8 

% 


% 

R 9 

0.0470 

-0.0208 

-0.0488 

0.0037 

1.094 

0.1464 


0.0376 

-0.1491 



Ring 

C - Forces in lb. 

Moments in 

in. 

-lb. 


-0.0297 

*-0.7775 

0.0960 

-0.0294 

-1.957 

-0.0864 


- 0.0501 

0.1454 

■ 


Ring 

D - Forces in lb. 

Moments in 

in. 

-lb. 


-0.0216 

-I .2798 

-0.0275 

-0.0252 

-1.7565 

JQ . 1682 


-0.0441 
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NACA TN No. 1014 Fig. I 

(a) ACTUAL CYLINDER 



STRINGER SECTION B-B 


(b) SIMPLIFIED CYLINDER 



RING PLUS EFFECTIVE SHEET 
SECTION C-C 


FIG. I. MONOCOQUE CYLINDER 

SHEET : 24 ST ALCLAD 
REINFORCEMENTS : 24 ST ALUM. ALLOY 




FIG. 2 UNIT AXIAL DISPLACEMENT 
OF A CORNER OF A PANEL 



FIG. 3 SIGN CONVENTIONS 


NACA TN No. 1014 


Figs. 4,5,6 




FIG. 4 UNIT TANGENTIAL DISPLACEMENT 


/ 



FIG. 5 UNIT RADIAL DISPLACEMENT 



FIG. 6 


UNIT ROTATION 





Figs. 7,8 


NACA TN No. 1014 


FORCES AND MOMENTS ACTING ON CONSTRAINTS. FOR SIGNS SEE FIG. 10a. 


t a“' ( *Ai 

Wd, 


X B = 


n b — (« n rQ r +(a n ru II 


*c = FI - 1 

Tc=KF) r 

N c ={« n ru r 


X 0 = ‘(EA/a) D pi\jj.-/L m Xp“ -(EA/a) or -{EAAJ FH fl^ 


-D 

R 0 =-(ix r n Ii: -(« r r) m 


"Fp“ - (<x^.r) I +(a (& ^r) m + Ca^r^ 


X H (EA/a) FV ffL I A-^ 

VKD r - (cc t r )jg- 

Ru**(cc ».r*)yC"(.OCr.V" l ') T vr 


R^— (oc^r’^u - (« r r)jjj. Rp” - Coc t .r)fC‘(cc ri r)jjHc< t . 1 r^jjj — (wprijju r^*(c< r r’)- t -^~(oc r p') T cr 

fcc ri F'-'i m Np* - (cCnPL'jjHa n rL') jj-{cc n r Lj ^ (a n f N H ^a n rUi - (« r ,r'-ilz 

h • i 


Xj“ Fl.ju * 

Rj“ (cc TV 
Nj-CanFOm 
* J 


*k“ Ajjf* - Ajy 

R^-OXfT^+CaprVsa; 

N K --(a^LyKa n rLV 


Xp* J*l-Txr 
T p — ( cc^Djj- 
Rp*~fcCrF^1V 


FIG. 7. EFFECT OF UNIT AXIAL DISPLACEMENT OF F. 
r=4t jl-^L 


FORCES AND MOMENTS ACTING ON CONSTRAINTS. FOR SIGNS SEE FIG. 10 a. 


V-KHn 

VK% 

^/C^n^-tAUix 


VK^VKAi 

V (tt'! BF -(a^ I -(ocyO II • 

Rb— (rUB^-Coc^oc^A^+CoCpOCtAlu. 

Ng- (^Bp-Coc^o^AL^-Coc^tADu- 


Xc’KD.i 

T c - («tA-\ 

N & -(a n °c t A.L\ 


X 0 'KrvK r V 


x F = ~WFV’K r W _ (o < tr )m + (oc^r*)jj- 
Tp^^^t^^Ctt^^Co^A^Co^A^^cc^A^jjjKoc^A^ 


Tf (oc|A’) I +(a^A'ljs r 


MvvMiE (k p «iA’W R H - (fVSAV(*r<*tAW 

Ntffe^lVi^ALV Np“^^ MB ^^ MF j^oc^^LjK« n o(^L^-(c( n ccAL\2^« n c( t Ai)j Z N H -(a n oc t AU I + -(oc n « t AU I 


X|^* (^n nx ”■ to v ^ 

Ty^ (tt) p K -(«^A'l jjj- (oc^Mjj 

R|C ^ VK~^r° < tAW"^r K tA^IZ 

N K = ^Vk- ( a n OC^A-Lljjf (eC^oCy ALijjj 


n,in ttn q i,ig: 

FIG. 8. EFFECT OF UNIT TANGENTIAL DISPLACEMENT OFT. 
F=-%L A=£ik 



Figs. 9,10 


NACA TN No. 1014 


FORCES AND MOMENTS ACTING ON CONSTRAINTS. FOR SIGNS SEE FIG. 10 a. 

A B C 




' * 


b 




FIG. ia RADIAL DEFLECTIONS FOR STRINGERS 
DUE TO PURE BENDING. 

M- 8155.4 IN. LB. 

POSITIVE SIGN INDICATES INWARD DEFLECTION 


FIG. II. AXIAL DEFLECTION OF RINGS 

PURE BENDING M - 6155.4 IN. LB. 
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NACA TN No. 1014 Figs. I 1,12 
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FIG. 15. AXIAL DEFLECTION OF RINGS 

COMPRESSION P - -1286 ' LB. 



FIG. 16. RADIAL DEFLECTIONS FOR STRINGERS 
DUE TO PURE COMPRESSION. 
R--I286 LB. 


POSITIVE SIGN INDICATES INWARD DEFLECTION 
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Figs. 15,16 NACA TN No. 1014 




DISTANCE FROM HORIZONTAL DIAMETER, INCHES 




nr 



35,000 IN. LB. 
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* V 




Figs. 19,20 NACA TN No. 1014 












DISTANCE FROM HORIZONTAL DIAMETER , INCHES 



FIG. 24. SHEAR NG STRESS DUE TO PURE BENDING 
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FIG. 25. SHEARNG STRESS DUE TO PURE COMPRESSION 


M *> 35,000 IN. LB. 


P » -1286 LB. 




r 




Figs. 24,25 ■ NACA.TN No. 1014 



MAXIMUM STRESS XKT 3 
LB. PER SQ. IN. 


NACA TN No. 1014 


Figs. 26,27 








Figs. 28,29 


NACA TN No. 1014 




FIG. 29. BENDING STRESS IN RING FOR SMALL CYLINDER. 

M = 42 000 IN. LB. 


1 





1 





FIG. 32 SOLUTIONS OF THE SIMPLIFIED CYLINDER 
COMPLETE FIELD 

M= 35,000 IN. LB. 



CUT FIELD 

M- 33,000 IN. LB. 



